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• Abstract 

We analyze lepton flavour violation (LFV) , as well as generation of the observed 
Oh. baryon-antibaryon asymmetry of the Universe (BAU) within a generalized mini- 

(— I I mal lepton flavour violation (MLFV) framework where we allow for CP violation 

both at low and high energies. The generation of BAU is obtained through ra- 



X 



diative resonant leptogenesis (RRL), where starting with three exactly degenerate 
right-handed neutrinos at Aqut, we demonstrate explicitly within the SM and the 
' MSSM that the splittings between their masses at the see-saw scale M^, generated 

by renormalization group effects, are sufficient for a successful leptogenesis for Miy 
even as low as 10^ GeV. The inclusion of flavour effects plays an important role 
in this result and can lead to the observed BAU even in the absence of CP vio- 
lation beyond the PMNS phases. The absence of a stringent lower bound on Miy 
in this type of leptogenesis allows to easily satisfy present and near future upper 
bounds on /i ^ 67 and other charged lepton flavour violating (LFV) processes 
even for Alfv = O(lTeV). We find, that the MLFV framework in the presence 
of heavy right-handed neutrinos and leptogenesis is not as predictive as MFV in 
the quark sector and point out that without a specific MLFV model, there is a 
rich spectrum of possibilities for charged LFV processes and for their correlation 
with low energy neutrino physics and the LHC physics, even if the constraint from 
the observed BAU is taken into account. While certain qualitative features of our 
analysis confirm findings of Cirigliano et al., at the quantitative level we find phe- 
nomenologically important differences. We explain the origin of these differences. 



1 Introduction 



One of the attractive and predictive frameworks for the description of flavour changing 
processes in the quark sector is the so-called Minimal Flavour Violation (MFV) hypothe- 
sis [T],[2] in which the Standard Model (SM) quark Yukawa couplings are the only sources 
of flavour changing and in particular CP-violating processesjj 

If only one Higgs doublet is involved in the spontaneous breaking of the underlying 
gauge symmetry, all flavour changing charged and neutral current processes are governed 
in the MFV framework by the CKM matrix [1] and the relevant local operators are 
only those present in the SM. As demonstrated in [5], the existing data on B^^ — B^^ 
mixing, Ek, B — > ^^7, B Xgl^l" and — s> tt^vv and the value of the angle (5 
in the unitarity triangle from the mixing induced CP asymmetry in i? ^ ipKs imply 
within this framework very stringent bounds on all rare K and B decay branching ratios. 
Consequently, substantial departures from the SM predictions are not expected if MFV 
with one Higgs doublet is the whole story. 

If two Higgs doublets, like in the MSSM, are involved and the ratio of the corresponding 
vacuum expectation values V2/V1 = tan/5 is large, significant departures from the SM 
predictions for certain decays are still possible within the MFV framework [2] in spite of 
the processes being governed solely by the CKM matrix. The most prominent examples 
are the decays Bd^s — A*'*'/^" with a subset of references given in [6] . The prime reason 
for these novel effects is the appearance of new scalar operators that are usually strongly 
suppressed within the SM and MFV models at low tan (3, but can become important and 
even dominant for large tan/3. The improved data on B^^g expected to come in 

this decade from Tevatron and LHC, will tell us whether MFV models with large tan (3 
are viable. 

One of the important virtues of the MFV in the quark sector are the relations [1], [7] 
between the ratios of various branching ratios and the CKM parameters measured in 
low energy processes that have universal character and are independent of the details of 
the specific MFV model. An example is the universal unitarity triangle common to all 
MFV models [1]. But also the fact that each branching ratio can be expressed in terms 
of the CKM parameters and quark masses measured at the electroweak scale or lower 
energy scales makes this scenario to be a very predictive framework. Moreover, neither 

^For earlier discussions of this hypothesis see [3]. 



fine tuning nor the introduction of unnaturally high scales of new physics are required 
to make this scenario consistent with the available data. 

The MFV scenario in the quark sector in question, although simple and elegant, suffers 
from the following problem. In the absence of new complex phases beyond the CKM 
phase, it cannot accommodate the observed size of the baryon asymmetry of the universe 
(BAU) to be denoted by 77^ in what follows. The question then arises, whether one could 
still explain the right size of rjB within the MFV context by considering simultaneously 
the lepton sector, where the BAU can in principle be explained with the help of leptoge- 
nesis [El [9]. While this is the most natural possibility, other directions could be explored 
in principle. 

Before addressing this question let us summarize what is known in the literature about 
the MFV in the lepton sector. Last year, Cirigliano, Grinstein, Isidori and Wise [10] 
in an interesting paper formulated MFV in the lepton sector (MLFV) both with the 
minimal field content and with the extended field content, where three degenerate right- 
handed heavy neutrinos u}^ with masses M* are added to the SM fields and the see-saw 
mechanism [TT] is responsible for the generation of the light neutrino masses with the 
see-saw scale denoted by M^, in what follows. Analyzing charged lepton flavour violating 
(LFV) processes, like /i — 67 and fi ^ e conversion in nuclei in these two scenarios in 
the absence of CP violation, they reached two interesting conclusions: 

• Measurable rates for LFV processes within MLFV are only obtained when the scale 
for total lepton number violation (Aln = C>{M,y)) is by many orders of magnitude, 
typically a factor 10^ — 10^, larger than the scale of charged lepton flavour violation 
(Alfv). 

• Similarly to MFV in the quark sector, the ratios of various LFV rates like ^ 
e'-f)/B{T fi'-f) are unambiguously determined in terms of neutrino masses and 
mixing angles measured in low energy processes. 

Various phenomenological aspects of MLFV, as formulated in [10], have been subse- 
quently discussed in [T2] . 

The MLFV framework in [10] does not include CP violation, neither at low energy nor 
at high energy, being a necessary ingredient in the generation of the BAU. Moreover, 
possible renormalization group effects between the low energy scale 0{Mz) and the high 
energy scales, like the see-saw scale M^, and the GUT scale Aqut, have not been taken 
into account in [10]. It is then natural to ask: 
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• whether a successful leptogenesis is at all possible within a MLFV framework in 
which flavour violation is governed solely by Yukawa couplings, 

• how the findings of [10] are modified, when CP violation at low and high energy 
and the renormalization group effects (RGE) in question are taken into account, 

• whether a successful leptogenesis in the MLFV framework puts stringent con- 
straints on charged LEV processes. 

The main goal of our paper is to answer these three questions. In fact, as we will 
demonstrate explicitely in Section 5, it is possible to obtain the correct size of r]B in 
the MLFV framework with three heavy right-handed neutrinos that are assumed to be 
degenerate in mass at Aqut- Other choices for this scale could be considered but Agut 
seems to be the most natural one. The breakdown of this degeneracy through RGE, 
that are governed by Yukawa couplings, combined with new sources of CP-violation in 
the heavy neutrino sector allows to obtain the correct size of 775 in the framework of the 
resonant leptogenesis in particular when flavour effects are taken into account. As this 
type of leptogenesis is generated here radiatively and not put by hand as done in most 
literature sofar, we will call this scheme radiative resonant leptogenesis (RRL) in what 
follows. 

The fact that within the MLFV framework one is naturally led to RRL, has significant 
implications on charged LFV processes, which could in principle be used to distinguish 
this scenario from other extensions of the SM. In particular, while other types of lepto- 
genesis, with hierarchical right-handed heavy neutrinos, imply generally rather stringent 
lower bounds for the lightest u}^ mass, in the ballpark of (9(10^ GeV) or higher, the 
values of M^y in RRL are allowed to be by many orders of magnitude lower. As the 
branching ratios for /j Ij'-f are proportional to M^/A^py [lO], it is relatively easy to 
satisfy the present and in the near future available upper bounds on these processes 
by simply choosing sufficiently small value of M^. Conversely, by choosing Mj^ to be 
larger than say 10^^ GeV, it is in principle possible to obtain the values of B{fi — cj) 
close to expected bounds from PSI even if Alfv is as high as 100 TeV. This means that 
non-observation of /i ^ 67 with the rate 10~^^ at PSI will not necessarily imply within 
the general MLFV framework that Alfv is very high. Conversely, the observation of 
/i — > 67 will not necessarily imply LFV physics at scales (9(1 TeV). In other words, 
without a specific MLFV model there is a rich spectrum of possibilities for charged LFV 
processes within the general MLFV framework, even if the constraint from r]B is taken 
into account. 

Thus one of the main messages of our paper is the realization that the MLFV framework 
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in the presence of heavy right-handed neutrinos and leptogenesis is clearly not as predic- 
tive as MFV in the quark sector. This is also related to the fact that new physics, even 
lepton conserving one, that could be present between energy scales Mz and M;^, could 
have in principle an important impact on various observables, like B{li Ijl)-, through 
RGE. 

In the advanced stages of our project a paper by Cirigliano, Isidori and Porretti [13] ap- 
peared, in which the idea of the incorporation of leptogenesis into the MLFV framework 
has been put forward in the literature for the first time and its implications for charged 
LFV processes have been analyzed in detail. While in agreement with the general pre- 
dictions of [13] we find that successful leptogenesis is possible for high scales My > 10^^ 
GeV, and contrary to that paper our detailed numerical analysis demonstrates that this 
is also true for much lower scales, weakening the implications for charged LFV processes 
found in that paper. Most importantly we do not confirm the lower bound of 10^^ GeV 
for My found by these authors which has significant implications for charged LFV pro- 
cesses as stressed above. The inclusion of flavour effects in the leptogenesis in our paper, 
that has been left out in [13] and the use of approximate formulae in that paper as op- 
posed to a full numerical analysis present here, brings in significant differences in these 
two analyses for My < 10^^ GeV. We will summarize the agreements and differences 
between [13] and us in Section [531 

At this stage it is worth also mentioning that there may be other equally reasonable 
definitions of MLFV. In this paper, we will only consider a conservative generalization 
of the initial proposal for MLFV [10]. However, it is clear that one may have other well 
motivated but different proposals for MLFV. In particular one should keep in mind that 
within the seesaw mechanism neutrinos acquire a mass in a manner which differs signif- 
icantly from the one in the quark sector. In fact it has been suggested [H] that the fact 
that neutrino masses arise from the seesaw mechanism is the key point in understanding 
why leptonic mixing is large, in contrast with small quark mixing. Therefore a reasonable 
definition of MLFV may differ from MFV in the quark sector. In our opinion, only in the 
presence of a theory of lepton flavour, where Yukawa couplings would be constrained by 
family symmetries, can one define in a unique way what MLFV should be. The question 
of different definitions of MLFV has been recently addressed in an interesting paper by 
Davidson and Palorini [T5] . 

Our analysis involves several points and it is useful to list them one by one already at 
this stage. 

• As already stated above, in the framework of the MLFV the right-handed heavy 
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neutrinos are assumed to be degenerate in mass at some high energy scale in order 
to exclude possible new sources of flavour violation. However, the exact degeneracy 
of M* is not RG invariant and can only be true at a single scale which we choose to 
be Agut = C(10^^ GeV). RGE between Agut and the see-saw scale <^ Agut 
break the degeneracy between M* at M^, a welcome result for leptogenesis that 
vanishes in the limit of degenerate M*. This structure is the basis of the so called 
radiative leptogenesis [Ml [17] that has been first considered in the case of two 
degenerate neutrinos in [T6 l fTT t ITS]. An important ingredient of this framework is 
the resonant leptogenesis [191 EOl EI] . Therefore we will call this framework RRL 
as stated above. Our analysis is one of the first that considers the case of three 
degenerate neutrinos and includes flavour effects in RRL. 

• As the values of the light neutrino masses and of the parameters of the PMNS 
mixing matrix [22], that enter the formulae for charged LFV processes, are not 
to be evaluated at the low energy scale but at the high energy scale M^, the 
MLFV relations between neutrino masses, mixing angles and rates for charged 
LFV processes presented in [10] can be in principle significantly modified through 
the RGE between the Mz and M^, scales, changing the conclusions about the value 
of the ratio Mj^/Alfv necessary to obtain visible charged LFV rates. While it is 
conceivable that in certain MLFV scenarios RGE could be neglected, the example 
of the MSSM with a large tan/3, presented in this context in [23], shows that 
the RGE in question could in principle modify — » Ijj) by a few orders of 
magnitude. 

• The requirement of a successful BAU with the help of leptogenesis and in fact in 
general, necessarily brings into play CP violation. Neglecting flavour effects in the 
Boltzmann equations, the relevant CP violation is encoded in a complex orthogonal 
matrix R in the parameterization of Y^, by Casas and Ibarra |24j . 

As analyzed already in several papers in the context of supersymmetric models, 
the size of the imaginary parts of R, crucial for generating the observed BAU in the 
framework of leptogenesis, can change the rates of charged LFV processes by several 
orders of magnitude. See, in particular [25], but also p3l [2il [26l [271 [28| [3T1 [M] . 
We note that when flavour effects are important, the generation of the BAU could 
be possible without complex phases in R [l6]. 

• The inclusion of CP violation at low energy with the help of the non-vanishing 
phase 5pMNS [22] has only a moderate impact on the results in [10] but in the 
presence of a complex matrix R (see above) non- vanishing Majorana phases in the 
PMNS matrix can modify the results for LFV processes in [10] both directly and 
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indirectly through RGE mentioned above. The numerical studies in [26l [271 123] 
show that such effects can be in principle significant. 

Our paper is organized as follows. In Section [2] we present the generalization of the 
formulation of MLFV given in [10] by including low energy CP violation in the leptonic 
sector with the help of the PMNS matrix [22] and the high energy CP violation necessary 
for the leptogenesis of BAU. The parametrization of the neutrino Yukawa coupling 
of Casas and Ibarra [24] turns out to be useful here. 

In Section [3] we analyze the issue of the breakdown of the mass degeneracy of heavy 
neutrinos by radiative corrections. This breakdown is necessary for leptogenesis to work 
even if R is complex. Assuming then this scale to be the grand unification (GUT) 
scale, we discuss the renormalization group equations in the SM and the MSSM used to 
generate the splitting of Mj at scales 0{M^), where the heavy neutrinos are integrated 
out. The results of this section are a very important ingredient of the leptogenesis that 
we consider in Section H] and in particular in [51 

In Section [H as a preparation for Section [5l we present some numerical aspects of the 
flavour changing radiative charged lepton decays k Ij'y and of the CP asymmetries in 
the right-handed neutrino decays. 

In Section [5l the most important section of our paper, we describe the scenario of radia- 
tive resonant leptogenesis in the case of three quasi-degenerate right-handed Majorana 
neutrinos. In this context we include in our analysis recently discussed flavour effects 
that definitely cannot be neglected. The main result of this paper is the demonstration 
that the right value of tjb can be obtained in this framework. A plot of tjb versus M^, 
demonstrates very clearly that already for M^, as low as 10^ GeV leptogenesis becomes 
effective and that flavour effects are important. We compare our results with existing 
literature and explain why in contrast to [13] we do not find a stringent lower bound on 
M,. 

Finally, we return to the k — > /j7 decays and use the knowledge collected in Sections [3l 
and [5l to present a brief numerical analysis of yU ^ 67 that illustrates the points made 
above. We restrict our analysis to tan/5 < 10 so that RGE between Mz and are 
small and other effects can be transparently seen. 

In Section [531 we compare our analysis and our results with [13]. We conclude in Sec- 
tion [6l 
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2 Basic Framework 



2.1 Preliminaries 

The discovery of neutrino oscillations provides evidence for non- vanishing neutrino masses 
and leptonic mixing, leading to lepton-flavour violation. In the SM, neutrinos are strictly 
massless, since Dirac masses cannot be constructed due to the absence of right-handed 
neutrinos, and left-handed Majorana masses are not generated due to exact {B — L) 
conservation. 

The simplest extension of the SM which allows for non-vanishing but naturally small 
neutrino masses, consists of the addition of right-handed neutrinos to the spectrum of 
the SM. This extension has the nice feature of establishing on the one hand a lepton quark 
symmetry and on the other hand being naturally embedded in a grand unified theory 
like 5*0(10). Since right-handed neutrinos are singlets under U{1) x SU{2) x S'?7(3), 
Majorana neutrino masses should be included, with a mass scale My which can be 
much larger than the scale v of the electroweak symmetry breaking. Apart from Mr, 
Dirac neutrino mass terms are generated through leptonic Yukawa couplings upon 
gauge symmetry breaking. The presence of these two neutrino mass terms leads, through 
the seesaw mechanism |TT], to three light neutrinos with masses of order v"^ /M^ and three 
heavy neutrinos with mass of order M^. The decay of these heavy neutrinos can play 
a crucial role in the creation of a baryon asymmetry of the universe (BAU) through 
the elegant mechanism of baryogenesis through leptogenesis P, E]. In the presence of 
neutrino masses and mixing, one has, in general, both CP violation at low energies which 
can be detected through neutrino oscillations and CP violation at high energies which is 
an essential ingredient of leptogenesis. The connection between these two manifestations 
of CP violation can be established in the framework of specific lepton flavour models. 

In this paper, we study lepton-fiavour violation in this extension of the SM, assuming 
minimal lepton flavour violation (MLFV) but allowing for CP violation both at low and 
high energies. The case of no leptonic CP violation either at low or high energies, was 
considered in [10] where the suggestion of MLFV was first presented. The first discussion 
of CP violation at low and high energy in a MLFV framework has been presented recently 
in [13]. We will compare the results of this paper with ours in Section [531 
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2.2 Yukawa Couplings and Majorana Mass Terms 



We add then three right-handed neutrinos to the spectrum of the SM and consider the 
following leptonic Yukawa couplings and right-handed Majorana mass terms: 

l^Y = -eRYE(t)^LL - i^rY^cPLl + h.c. (1) 

= -^I^rMrUr + h.c. , (2) 

where Ye, Y^ and Mr are 3x3 matrices in the lepton flavour space. In the limit 
£y = Cm = the Lagrangian of this minimal extension of the SM has a large flavour 
symmetry 

SU{?>)l X SU{?>)e X SU{?>),^ X U{1)l x U{1)e x f/(l).^, (3) 

which reflects the fact that gauge interactions treat all flavours on equal footing. This 
large global symmetry is broken by the Yukawa couplings Ye, Yy and by the Majorana 
mass terms Mr. A transformation of the lepton fields: 

Ll -> VlLl, cr VeCr, ur Vy^UR (4) 

leaves the full Lagrangian invariant, provided the Yukawa couplings and the Majorana 
mass terms transform as: 

n ^ r; = Vy.YyVl (5) 
Ye^Ye = VeYeVI (6) 
Mr M'^ = V:^MrVJ^, (7) 

which means that there is a large equivalent class of Yukawa coupling matrices and 
Majorana mass terms, related through (jSD-dZD, which have the same physical content. 
The MLFV proposal [10] consists of the assumption that the physics which generates 
lepton number violation, leading to Mr, is lepton flavour blind, thus leading to an exactly 
degenerate eigenvalue spectrum for Mr, at a high-energy scale. As a result, in the MLFV 
framework, the Majorana mass terms break S'f/(3);^^ into 0(3),^^. 



2.3 Leptonic Masses, Mixing and CP Violations 

Without loss of generality, one can choose a basis for the leptonic fields, where Ye and 
Mr are diagonal and real. In this basis, the neutrino Dirac mass matrix niE, = vY^ is an 
arbitrary complex matrix, therefore with nine moduli and nine phases. Three of these 
phases can be eliminated by a rephasing of Ll. One is then left with six CP violating 
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phases. There are various classes of phenomena which depend on different combinations 
mu, mj), m|) or equivalently Y^, Yj and Yj: 

A) Leptonic mixing and CP violation at lov^^ energies: Since we are working 
in the basis where the charged lepton mass matrix is diagonal and real, leptonic 
mixing and CP violation at low energies are controlled by the PMNS matrix U„ 
[22], which diagonalizes the effective low energy neutrino mass matrix: 

U'^{m^)effUy = dy, (8) 

where = diag(mi, m2, ms), with rrii being the masses of the light neutrinos and 

m 

{m,),ff = -v^YjD-^'Y,, (9) 

where denotes the diagonal matrix Mr and v = 174 GeV. In the case of MLFV, 
Dfj = Myl and one obtains at My ^ Aln 

K)e// = -^/^^Yy. (10) 

Consequently YjY^ is the quantity that matters here. 

B) Lepton flavour violation: The charged LFV depends on the other hand on the 
combination YjY^ with Yy again normalized at the high energy scale M^. We will 
return to this point in Section HI 

C) CP violation relevant for leptogenesis: The generation of BAU through lep- 
togenesis starts by the production of a lepton asymmetry which is proportional to 
the CP asymmetry in the decays of heavy Majorana neutrinos. This CP asym- 
metry involves the interference between the tree-level amplitude and the one-loop 
vertex and self-energy contributions. It has been shown [30] that the CP asym- 
metry depends on the neutrino Yukawa couplings through the combination Y^Yj . 
Again as in classes A and B, Y",^ is evaluated here at the scale M^. When flavour 
effects in the Boltzmann equations become important, the non-summed products 
(Yy)ik{Yi,)*f, corresponding to different lepton flavours k can attain relevance. 

2.4 An Useful Parametrization 

In order to analyze in a systematic way the above phenomena and study the implied rela- 
tions among low-energy lepton mixing data, lepton flavour violation and leptogenesis in 
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different scenarios classified below, it is convenient to choose an appropriate parametriza- 
tion for Y^. We use the following parametrization [23] of the neutrino Yukawa couplings: 

(/D^)-iy; = ^i?v^f/t, (11) 

where R is an orthogonal complex matrix {R^R = RR^ = 1), di, = diag(mi, m2, ma) 
and Dr = diag(Mi, M2, M3). 

It is instructive to count next the number of independent parameters on both sides of 
(fTTl) . The left-hand side of ( fTTl) is an arbitrary 3x3 complex matrix with nine real 
parameters and six phases, since three of the initial nine phases can be removed by 
rephasing Ll- It is clear that the right-hand side of ffTTl) also has nine real parameters 
and six phases. Indeed, R, dy and Uy have each three real parameters and moreover R 
and Uy have in addition each three phases. We consider now the case where the right- 
handed neutrinos are exactly degenerate, i.e. = M^l. We will show that three of 
the real parameters of R can be rotated away. Note that any complex orthogonal matrix 
can be parametrized as 

R = e^^e'^\ (12) 

with A12 real and skew symmetric. Now in the degenerate case an orthogonal rotation 
of z/R — > OrVr leaves the Majorana mass proportional to the unit matrix and defines 
a physically equivalent reparametrization of the fields ur. Choosing Or = e^^ we see 
immediately that 

n - O^rY, = ^ e-^^ R^uUl = ^ e*^^ v^[/t , (13) 
which shows that the physically relevant parameterization is given by Rdeg = e^^^- 

Using the parameterization in f|TT]) one finds that the matrix YjY^, which controls low- 
energy CP-Violation and mixing can be written as follows 

YjY, = -\{Ulf ^,R^ DrR^^UI = -^{UlfdyUl (14) 

V f 

where in the last step we have set Dr = M^l. 

On the other hand, the matrix YjY^y which controls charged LFV, can be written as 
follows (see also 1311) 



YjYy = ^ ^/d', R^ DrR ^yUl = ^ U, ^/lyR^R^.Ul . (15) 



Finally, the matrix Y^}, which enters in leptogenesis when flavor effects are not relevant 
is given by (see also [3T]): 

Y,Yl = \^RRd,R^ ^R = ^Rd,Ry (16) 
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We note that Y^Y^, depends only on f/y and d^, while Yj^Yj relevant for the leptogenesis 
only on d^, and R. This means that CP violation at low energy originating in the complex 
Ui, and the CP violation relevant for leptogenesis are then decoupled from each other 
and only the mass spectrum of light neutrinos summarized by d,y enters both phenomena 
in a universal way. 

In this respect the charged LFV, represented by f|T5|) . appears also interesting as it 
depends on d^, and R and consequently can also provide an indirect link between low 
energy and high energy CP violations and generally a link between low and high energy 
phenomena. 

2.5 Classification 

Having the parametrization of Y^ in ( JTlll at hand we can now spell the difference between 
the analysis of [10] and ours in explicit terms. Indeed, from the above considerations, 
it follows that possible relations among phenomena A,B,C, discussed in Section 2.3, 
crucially depend on the assumptions one makes about leptonic CP violation at low 
energies, as well as at high energies. One may consider then separately the following 
four scenarios: 

• Case 1: No leptonic CP violation either at low or high energies. The limit that 
all complex phases vanish leads to 

-YJY, = YX- (17) 

This is the case considered in [10], where a close connection is obtained between 
experimental low energy data on lepton mixing and the pattern of various charged 
LFV processes, that is the correlation between phenomena A and B in the absence 
of CP violation. It corresponds to choosing R and Up real. However, even in this 
case the RGE between the low energy scale at which the light neutrino masses and 
mixings are measured and the scale My at which Yy is evaluated could have an 
impact on the correlation in question. 

• Case 2: Leptonic CP violation at low energies, but no CP violation relevant 
for leptogenesis (barring flavour effects). This corresponds to assuming that the 
leptonic mixing matrix Uy contains CP violating phases so that YyYy is complex, 
but YyY^ is real or equivalently as seen in ( fT6|) R is real. 

• Case 3: CP violation relevant for leptogenesis but no low energy leptonic CP 
violation. This corresponds to having YyY} and R complex, but Uy real. 
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• Case 4: There is leptonic CP violation both at low and high energies, that is both 
f/j, and R are complex quantities. It should be stressed that Case 4 is of course 
the general case and, in fact, the most "natural" one, since once CP is violated by 
the Lagrangian the six CP violating phases contained in lead in general to CP 
violation both at low and high energies. 

2.6 Final Remarks 

It is clear that f|T5l) . depending on f/j,, R, and enables one to analyze separately 
the four cases considered here. In each case there will be simultaneously implications 
for lepton flavour violations, leptogenesis and low energy CP violation and mixing with 
certain correlations between them. These correlations can be affected by RGE between 
the low energy scale and M^. 

At this stage the following comments are in order: 

• Uy is relatively well known from oscillation experiments with the exception of S13, 
the phase 6 and the Majorana phases a and (3. In order to use it for the calculation 
of Yjj it has to be evolved by RG equations to My. 

• With the measured two mass differences squared from solar and atmospheric os- 
cillation data, the diagonal matrix is a function of a single parameter that we 
choose to be the mass of the lightest neutrino. Again these parameters have to be 
evaluated at the scale with the help of renormalization group techniques. 

• The matrix R depends on three complex parameters that influence simultaneously 
lepton flavour violation and leptogenesis as seen in (ITSll and (ITGll . respectively. 
Some constraints on R can then be obtained from these two phenomena but a 
complete determination of this matrix is only possible in an underlying theory 
represented usually by special texture zeros of Yy. 

• Finally, My can be restricted from the BAU in the context of the seesaw mechanism 
and if the eigenvalues of the right-handed neutrino matrix are hierarchical, the 
absolute lower bound on the lowest Mi is (9(10^) or even higher. In the case of 
the MLFV considered here the right-handed heavy neutrinos have to be quasi- 
degenerate in order to avoid new flavour violating interactions. In this case BAU 
can be explained with the help of RRL which combines the resonant leptogenesis 
considered in [201 EI] and radiative leptogenesis [161 lEl [E] • The lower bound on 
My can be significantly lowered in this case, as we will see explicitely below. 
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3 Radiative corrections in MLFV 



3.1 Preliminaries 

Our MLFV scenario defined in the previous section contains no free parameters beyond 
tlie neutrino masses, the PMNS matrix, a matrix of form Rdeg, an initial, universal heavy 
Majorana neutrino mass, and perhaps additional flavour-blind parameters that depend 
on the MLFV model. The rates for charged lepton flavour violation thus follow upon 
computing radiative corrections due to the degrees of freedom between the scales Mz 
and Aqut, and with suitable washout factors also the baryon asymmetry 77^. 

In this section we investigate how the CP- and flavour-violating quantities relevant to 
leptogenesis and charged lepton flavour violation, respectively, are radiatively generated. 
Since leptogenesis in the present framework can be considered as a generalization of the 
setup with two heavy singlets in [18] to the case of three degenerate flavours, we will also 
clarify what novelties arise in this case. This will be important in comparing our results 
to the existing literature. 

An important point will be that, due to the hierarchy between the GUT/flavour-breaking 
scale Agut and the neutrino mass scale Mi^, large logarithms appear such that the 
parameter counting for the coefficients Cj of flavour structures that has been recently 
presented in |T3] should be modified. Rather than being independent, the coefficients of 
structures containing different powers of Yukawa matrices are related by the renormal- 
ization group, while any additional independent effects are suppressed. Although this 
fact in principle increases the predictivity of MLFV, in our phenomenological sections 
it will still turn out insufficient to have correlations between high-scale and weak-scale 
observables. 



3.2 MLFV with a degeneracy scale 

We have defined our MLFV scenario to have a scale at which the masses of the right- 
handed neutrinos are exactly degenerate, such that the matrix Mr has no flavour struc- 
ture at all. In general, there will be additional flavoured particles in the theory. As a 
specific example, we consider the MSSM. Here the iVj are accompanied by heavy sneu- 
trinos N^, and there are also SU(2) doublet sleptons /j, transforming as 

l^Vj, N-^V*N' (18) 
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under the transformation (j4]). The Lagrangian then contains soft SUSY breaking terms 

= -Nrm%N^^ - i:m%l, + . . . , (19) 

where the elhpsis denotes further scalar mass matrices and trihnear scalar interactions. 
In general all matrices in Csoft have non-minimal flavour structure. The simplest general- 
ization of our degenerate scenario is then to extend the requirement of exact degeneracy 
to all mass matrices, similar to minimal supergravity. To be specific, we require all scalar 
masses to have the same value ttiq at the high scale and also require the A-terms to have 
the mSUGRA form A = aY with Y the corresponding Yukawa matrix and a a universal, 
real parameter of the theory. This example also provides us with a concrete value for 
the scale Alfv: LFV processes such as k — > Ij^y are mediated by loop diagrams involving 
sleptons and higgsinos or (weak) gauginos, and unless gaugino masses are very large, the 
scalar particles such as /j decouple at a scale A ~ mo. Hence the operators governing 
charged LFV are suppressed by powers of mo = Alfv- As in the case of the heavy 
Majorana masses, the generalized degeneracy requirement is not stable under radiative 
corrections, and for the same reason it is not renormalization scheme independent. 

3.3 Radiatively generated flavour structure and large logarithms 

As will be discussed in detail in the following section, the CP asymmetries necessary for 
leptogenesis require mass splittings between the decaying particles. The decaying parti- 
cles are on their mass shelo, but the degenerate initial conditions are usually specified 
in a massless schem^ (MS to be definite [32]). 

At one loop, the two mass definitions are related by a formula of the structure 

j^os ^ Mf^{fi) + CiMf^{fi) In — + nonlogarithmic corrections, (20) 

where fi ~ Aqut is the MS renormalization scale, q = 2(Y^Yj)ii/ {Wn'^) in the standard- 
model seesaw, and the nonlogarithmic corrections depend on our choice of massless (or 
any other) renormalization scheme. The resulting scheme dependence cannot be present 
in physical observables such as the BAU. Since this issue is usually not discussed in the 
literature on lepton flavour violation, let us elaborate on how it may be resolved. 

First, notice that while the nonlogarithmic terms in fl2Ul) are scheme dependent, the loga- 
rithmic corrections proportional to q are actually scheme independent. If IhAqut/Mu ^ 

^ We follow the treatment of [501 [H] (see also [12]), where sometimes the on-shell masses are replaced 
by thermal masses. (We will employ zero-temperature masses.) 

■^This is likely appropriate if the degeneracy is true to some flavour symmetry of an underlying theory, 
relating high-energy Lagrangian parameters and broken at the scale Aqut- 
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1, the logarithmic terms must be considered 0{1) and summed to all orders. This is 
achieved in practice by solving renormalization group equations. Similar resummations 
must be performed for all other parameters in the theory (such as Yukawa couplings). 
Correspondingly, the dominant higher-loop corrections to LFV observables and leptoge- 
nesis are approximated by using leading-order expressions with one-loop RGE-improved 
Yukawa couplings and masses. This is the leading-logarithmic approximation. Nonlog- 
arithmic corrections such as those indicated in ( l20l) are then sub-leading and should be 
dropped. 

What happens when the logarithms are not large is the following. If the MLFV frame- 
work is an effective theory for some fundamental theory where the degeneracy is enforced 
by a flavour symmetry, for instance the group (l3l), then the degeneracy holds in any 
scheme (that respects the symmetry) in the full theory and the scheme dependence ob- 
served in ( l20l) must be due to unknown threshold corrections in matching the underlying 
and effective theories. Since the flavour symmetry in MLFV, by definition, is broken 
precisely by the Yukawa matrices, this matching introduces all possible terms that are 
invariant under transformations (Hf5|6|7l) . A list of such structures has recently been 
given in [13], for instance. 



Mr = M, 



1 + c,{Y,Yj + (Y.Yjf) + C2{Y,YjY,Yj + (Y^YXyJf) + ■■■ ■ (21) 



The coefficients ci and C2 have been claimed by these authors to be independent 0{1) 
coefficients. Indeed these terms contain only non-logarithmic terms and (small) decou- 
pling logs when Mji is taken in the MS scheme, renormalized near the GUT (matching) 
scale. 

However, when computing the (physically relevant) on-shell in the case of Aqut 3> 
M^, large logarithms dominate both ci and C2. The leading logarithmic contributions 
are not independent, but are related by the renormalization group. C2 is quadratic in 
L = IylAcvt/M^, while ci is linear, and the RGE for Mr implies C2|l2 = |[ci|i]^. These 
logs are summed by RG-evolving to a scale fi ~ M,y. The additional conversion to 
on-shell masses is then again a sub-leading correction. 

Finally, we note that if there is no underlying symmetry, the degeneracy condition can 
again be true at most for special choices of scheme/scale, and must be fine-tuned. 

Numerically, the logarithms dominate already for mild hierarchies Aqut/^;/ > 10^, as 
then 2 In Aqut /M^ ^ 10. Let us now restrict ourselves to hierarchies of at least two 
orders of magnitude and work consistently in the leading-logarithmic approximation. As 
explained above, in this case non-logarithmic corrections both of the threshold type (in 
the coefficients q in ( !2T|) and in physical quantities (on-shell masses, CP asymmetries. 
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etc.) are sub- leading and should be dropped. In this regard our apparently "special" 
framework of initially degenerate heavy neutrinos turns out to be the correct choice at 
leading-logarithmic order. 

Finally we recall that the positions of the poles of the Ni two-point functions contain 
an imaginary part related to the widths of these particles. While not logarithmically 
enhanced, these are also scheme-independent at one loop (as the widths are physical), 
and it is unambiguous to include them in applications. In fact, these widths effects are 
often numerically important for the CP asymmetries in Ni decay [201 El], and we will 
keep them in our numerical analysis. 



3.4 Renormalization-group evolution: high scales 

For the running above the seesaw scale the relevant renormalization-group equations have 
been given in in [51] (in particular, last paper) for the SM and MSSM seesaw models. 
As the physical quantities studied below, such as leptonic CP asymmetries, involve mass 
eigenstates, it is convenient to keep the singlet mass matrix diagonal during evolution 
(see, e.g.. Appendix B of [52]): 

MRifx) = diag(Mi(/i),M2(/i),M3(/i)). 

Defining 

H = Y,Yl (22) 

and 



1 

167r' 

one obtains for the mass eigenvalues in the SM with right handed neutrinos: 



^ = T7^1iiW^gut), (23) 



^ = 2HuM, (no sum). (24) 

Note that due to the positivity of the right-hand side of ( 12^ . the running will always 
decrease the masses when running from the GUT to the seesaw scale. 

The matrix H satisfies the RGEs 

^ = [T,H]+3H'' -3 Y.Y^YeYJ + 2aH (SM), (25) 

^ = [T,H]+QH^ + 2 Y.Y^YeY^ + 2aH (MSSM) , (26) 
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where 

a 

3 



Tr(Fjn)+Tr(F;Fe) + 3rr(F„tyj + 3Tr(rir,)- (SM), (27) 



a = Tr{Y}Y,) + ^Tr{YlY^)--gl-^gl (MSSM), (28) 



ReH,, - z^lmH,, ^ j, SM), 



Mj-Mi " ^^-■■-■'t] " Mj +Mi 

T^'J = { -2 ^Rei/,, - 2 z ^M,, (z ^ j, MSSM), (29) 

(^ = j), 

and GUT normalization has been employed for gi. The matrix T satisfies t/ = Tf/, 
where M^\ij,) = U {nY M R{jj)U {jj) and m}j°^ satisfies the unconstrained RGEs given 
in [51]. Note that a is real and has trivial fiavour structure. Note the different relative 
signs in (l25l) and (l26l) : we will return to this point below. 



We now turn to a qualitative analysis of these equations and their impact on leptogenesis 
and flavour violation. Ignoring flavour effects in the Boltzmann evolution of charged 
leptons, the baryon asymmetry rjB is approximately proportional to the combinations 
Im((ifjj)^) = 2 ReHij ImHij {i ^ j), evaluated in the mass eigenbasis. At the scale Aqut, 
degeneracy of allows the use of an 5*0(3) transformation to make the off-diagonal 
elements of Rei7 vanish]^ As explained above, we should RG-evolve all parameters to the 
scale /i ~ Ml, to avoid large logarithms. Let us flrst consider the formal limit of vanishing 
charged lepton Yukawa couplings Ye for the SM case. It is instructive to split (|25l) into 
real and imaginary parts. The former satisfles 
dReH 



dt 



[ReT, ReH] - [ImT, \mH] + ?,^{ReHf - {\mHf^ + 2aReH. (30) 



To investigate how a nondiagonal ReH can be generated radiatively, assume that it is 
zero at some scale (initial or lower). Then ( l30ll reduces to 



dt 

(At t = 0, an extra term proportional to the offdiagonal part of (Imif )^ appears on the 
right-hand side of fl3Tl) .) Now evaluate this for the (2, 1) element and notice that Tij = 
and ImHij = for i = j. If there were only two heavy singlets in the theory, each term 
in each matrix product would require one (2, 1) element and one (1, 1) or (2, 2) element 
from the two matrix factors. For example, 

(ImT lmH)2i = ImTzi ImHu + ImTaa Imi/ai = 0, (32) 



To see this, notice that H is hermitian, so HeH is real symmetric. That is, it can be diagonahzed 
by a real orthogonal (and hence unitary) transformation of the right-handed neutrinos. Now if all three 
neutrinos are degenerate, such a rotation affects no term in the Lagrangian besides Y,y. 
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and similarly for the other terms. Consequently, 

Rei72i = ^ — = 0. (33) 

We see that there is no radiative leptogenesis in the two-flavour case when Ye = 0. This 
is consistent with the approximate equation (12) in [18], where Reif2i was found to be 
proportional to y'^. It is easy to see that the argument breaks down in the three- flavour 
case. For instance, 

{{lmHY)2i = lraH2ilmHn + Im^r22lmi/2i + Imi/sslmi/si = Imi/aslmi/ai, (34) 

which is in general not zero. The other terms in (!3T|) are also proportional to Imif23lni-^3i- 
We see that three generations of heavy neutrinos are necessary and sufficient to generate 
leptogenesis without help from charged lepton Yukawas. 

Once we restore the charged lepton Yukawas, they will also contribute. The impor- 
tant qualitative difference is that, whereas the contribution involving the charged- lepton 
Yukawas is only logarithmically dependent of the seesaw scale (as seen in eqs. fl58l) - (l60l) 
below for the two-flavour case, or from [21] for the three- flavour case), the pure Y^, contri- 
bution to the radiatively generated KeHij scales with because it contains two extra 
powers of Y^, as observed in the three flavour scenario studied in [T3] . 

In summary, we expect the following qualitative behavior for the BAU as a function of 

• For small Yi: (small M^), the dominant contribution to KeHij and hence to tjb 
should be due to Ye- r]B turns out to be weakly dependent on M^. 

• For large (large M^), in the three- flavour case there is a relevant contribution 
proportional to ((Imif)^)jj. Since it contains two extra powers of Y^, with respect 
to the contribution proportional to y"^, rj^ scales linearly with M^. 

• In the case of only two heavy flavours, rjB is weakly dependent on M^, over the 
whole range of M^. We will therefore include an "effective" two-flavour scenario in 
our numerical analysis. 

Let us stress that we reached these qualitative conclusions only upon neglecting flavour 
effects in the Boltzmann evolution of the products of the Ni decays. We will return to 
these points in Section H] and in Section [5l where we perform a detailed quantitative 
analysis. 

Finally, let us briefly discuss k Ij'-f. In MLFV these radiative lepton decays are 
governed by Aij = Y^Yy (and structures involving more powers of Yukawa matrices). 
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In the case of the SM, the rates are known to be essentially zero due to a near perfect 
GIM cancellation among the tiny neutrino masses. From the point of view of MLFV, 
this smallness can be traced to the fact that, in the SM, the LEV scale is equal to the 
LNV scale ~ M^. 

On the other hand, in the more generic case of the MSSM, there are additional contribu- 
tions mediated by slepton-higgsino or slepton-gaugino loops suppressed only by a scale 
Alfv ~ f^ii of order TeV, as discussed in Section 13.21 Linearizing the RG evolution, the 
charged slepton soft mass matrix acquires the form [50] 

m2(M,) = mil - L-i^iQml + 2al) + . . . , (35) 

where the dots denote terms governed by charged lepton Yukawa couplings Ye or con- 
serving lepton flavour. Note that the flavour structure in the soft terms is generated 
at a high scale and that, unlike the case of CP asymmetries in Ni decay, the necessary 
flavour structure A is already present at the initial scale Agut- Hence the RGE running 
of A merely gives a correction. Note also that there is dependence on the MLFV model 
beyond the choice of LEV scale due to the (in general unknown) RGE coefficients in (the 
relevant analog of) (!35|l . 

3.5 RGE evolution below M„: PMNS matrix and Aij 

So far we have ignored renormalization effects in equations such as ffTTl) . identifying d^, 
and Ui, with the physical (light) neutrino masses and mixing matrix, while the objects 
Yi, and Dji are defined at a high scale. However, to be orthogonal the matrix R has to 
be defined with all objects given at the same scale. Now it is well known that using low- 
energy inputs in can be a bad approximation because there are significant radiative 
corrections between the weak and GUT scales. However, as investigated in [52], both in 
the SM and in the MSSM with small tan j3 the main effect below Mi, is an approximately 
universal rescaling of the light neutrino masses. This results in larger magnitudes of the 
elements of Y^, extracted by means of ffTTj) but in a weak running of the matrix U,^. 
Above the scale M^,, even though the heavy singlets are now dynamical, one can still 
define an effective neutrino mass matrix through the seesaw relation Q). However, the 
evolution becomes more involved, as in the presence of heavy singlets there are additional 
contributions to the running involving Y^. To deal with this situation, where some of 
our inputs are specified at the weak scale, while the matrix Rdcg is defined at the scale 
Agut, we employ an iterative procedure detailed in Appendix |Al As was the case for the 
evolution above M^, also the RGE effects below M^, and consequently the relation of e.g 
YjYy to the input parameters necessarily depends on the details of the MLEV model. 
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4 Numerical Analysis: B(li 
metries in i^r decay 



Ij^) and CP asym- 



4.1 Preliminaries 



For our numerical analysis we take our input parameters at the weak scale, except for 
the matrix Rdcg, which has to be defined at the scale Aqut- From these inputs we find a 
consistent set of parameters at the seesaw scale M^, where the CP asymmetries as well 
as B{li — > Ij'-f) are calculated, through the iterative procedure given in Appendix [A] For 
the running we use the package REAP [53] . 
For the PMNS matrix we use the convention: 



( C12C13 
-S12C23 

\ S12S23 - Ci2C23Si3e 



■5l2Cl3 



Sl3e" 



Cl2S23Sl3e 
j<5 



if, 



iS 



C12C23 - Si2S23Sl3e S23C13 
-S23C12 - Si2C23Sl3e*'^ C23C13 J 



■ V 



(36) 



and V = Diag(e*"/^, e*^/^, 1) where a and P denote the Majorana phases and S denotes 
the Dirac phase. We parameterize the complex orthogonal matrix R as follows: 



(37) 





/ 


C12 


S12 


0\ 




/I 


\ 


/ Ci3 





Sl3^ 


R = 




-S12 


C12 










C23 S23 




1 







\ 














-S23 C23 / 


\-Si3 





Cl3 / 



with Sij = sin 6ij, with 6ij in general complex: 



Oij 



(38) 



In the degenerate case, the angles Xij can be made to vanish by a redefinition of the 
right-handed neutrinos, i.e. a matrix of the form R^eg is parameterized by three real 
numbers i/ij. 

In the following, we use maximal atmospheric mixing C23 = S23 = 1/ and a solar 
mixing angle 6*801 = 33°, with corresponding values for its sine s = S12 and cosine c = c^. 
For the sine of the CHOOZ angle S13 and the phases we allow the ranges 



< si3 < 0.25, < a, /3, 5 < 27r, 
and for the light neutrinos we use the low energy values 

- m? = 8.0 ■ 10"^ eV^ 



^"^sol 



Am 



atm 



jmg — m 



2.5 ■ 10"^ eV 



(39) 

(40) 
(41) 
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< m|^shte«t < 0.2 eV (42) 

with TTijjg'^t'^*^* = rriilm^) for normal (inverted) hierarchy, respectively. See [33] for a 
detailed discussion of the neutrino masses and mixing. For the heavy neutrino mass 
scale, we consider a wide range 

10^ GeV <M^< 10^^ GeV, (43) 

and the CP violating parameters yij are all taken in the range [—1,1] if not otherwise 
stated. 



4.2 Pert ur bat ivity bounds 

In the MLFV framework the magnitudes of the Yukawa couplings Y^, are very sensitive 
to the choice of M^, rri]^^^^'^^^ and the angles in the matrix Rdeg, as is evident from (fTTjl . 
To render the framework perturbative, we impose the constraint 

Att 

where y^i^^ is the largest eigenvalue of YjY,y . By means of (|T5l) . it translates into a bound 
on R)R = and the angles Uij that scales with and hence is most severe for a 
large lepton-number-violating scale. Analogous bounds apply to other dimensionless 
couplings whose number depends on the precise MLFV model. For instance, in the SM 
there is also the Higgs self coupling Xh-, whereas in the MSSM there is no such additional 
coupling. 



4.3 Lepton Flavour Violation and li Ij^ 

Following Cirigliano et al. [ID] we consider the normalized branching fractions defined 

as 

where B{li Ijj) is the true branching ratio and r^e = 1-0, = 5.61 and Tt-^ = 5.76. 
Assuming first the heavy right-handed neutrinos to be degenerate but not making the 
assumptions of i? = 1 and Ui, being real as done in [10] , the straightforward generalization 
of (29) in [TO] is 

B{h ^ /,7) = 3847rV-^|A,,f (46) 
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Here v = 174 GeV is the vacuum expectation value of the SM Higgs doublet |f| Alfv is 
the scale of charged lepton flavour violation, and C summarizes the Wilson coefficients 
of the relevant operators that can be calculated in a given specific model. They are 
naturally of 0{1) but can be different in different MLFV models. As we would like to 
keep our presentation as simple as possible, we will set \C\ = 1 in what follows, bearing 
in mind that in certain scenarios C may differ significantly from unity. Thus the true 
B{li — s> Ij'-f) can be different from our estimate in a given MLFV model, but as C is, 
within MLFV, independent of external lepton flavours, the ratios of branching ratios 
take a very simple form 



The most important objects in (146|) and (147|) are 



mn 



12 ■ 



(47) 



A,, = {Y}Y,),, = ^{U,^/d',R^DnR^Mlh, (48) 



which in the limit oi R = 1, Dr = M^l, and being real reduce to Ajj as given in 

(14) of [in]. 

With the formula PHI) at hand we can generalize the expressions for Ajj in (24) of 
[To] to the general case of complex R and U,^. To this end we will use the standard 
parametrization of the PMNS matrix Up in ( l36l) and the parametrization of R in ( !37|) . 
As the general expressions for Ajj in terms of Xij and yij are very complicated, we give 
in Appendix B explicit formulae setting all Xij = yij = except for yi2 ^ 0. We will 
discuss in our numerical analysis also the cases for which ?/i3 and 1/23 are non-vanishing. 
As mentioned above, setting Xij = is in accord with the degeneracy of the right-handed 
neutrinos. Once this degeneracy is broken by RG effects, the xij become non-zero. 

Recall from Section [3] that Ajj evolves above the scale My and the flavour structures it 
affects, such as the slepton mass matrix m?, also evolve between My and Alfv (and the 
resulting effective operators below Alfv also evolve). Moreover, the flavour- violating 
piece in, for example, m? is not exactly proportional to A at the scale My beyond 
leading order because these objects satisfy different RGBs between My and Agut- All 
this running depends, beyond the operator, also on the details of the model. Below 
the seesaw scale the flavour-non-universal contributions are governed by (although 
trilinear couplings such as the A-terms in the MSSM can also contribute), which is 
analogous to the case of the PMNS matrix. Based on the experience that the running of 



= \/vl + V2 for two-Higgs-doublet models such as the MSSM. Powers of sin/3 can be absorbed 
into C or into a redefinition Alfv ^j%v 
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Figure 1: Double ratios of k Ij'j for the MSSM with tan/3 = 2. Left plot: All 
parameters varied, right plot: no phases and only yi2 ^ 0. For a discussion, see the text. 



the PMNS angles is weak in the SM and the MSSM unless tan (3 (and hence t/r) is large, 
we ignore all these details and evaluate Aij at the scale M^. 

That Aij has to be evaluated at the high energy scale M^, and hence Ui, and di, have 
to be evaluated at by means of renormalization group equations with the initial 
conditions given by their values at Mz, has recently been stressed in particular in |31j . 
The dominant contributions to the flavour- violating pieces in the charged slepton masses 
matrix in the MSSM that is relevant for /j — >■ Ij'j are proportional to YjYi, and come 
from scales above M^, as seen for instance in equation (30) of [M] (where charged lepton 
Yukawas and A-terms have been dropped and only contribute at higher orders) and the 
fact that right-handed neutrinos and their Yukawa couplings are absent below that scale. 

All other parameters of a given MLFV model, hidden in the Wilson coefficient C in 
fH6|) . like slepton and chargino masses in the MSSM, would have to be evaluated at the 
electroweak scale and lower scales if a concrete value for C was desired. 

The ratio B{fi e'-f)/B{T fi'-f) is shown for the case of the MSSM with tan/5 = 2 
in Fig. [1] (left). All other parameters are varied in the ranges given above. We see that 
this ratio varies over about six orders of magnitude and B{fi — > 67) can be a factor 10'^ 
larger than B{t fi'y) in qualitative agreement with [21] • We have checked that 
the leptogenesis constraint, as discussed in Section [5l has no significant impact. This 
contradicts the findings of [13]. Even when constraining the Dirac and Majorana phases 
in the PMNS matrix to zero and allowing only for a single non-vanishing angle 1/12 at 
the scale Aqut, we can still have B{fi — > 67) ^ B{t fi^y). This is again in agreement 
with [2Sl[2n]- We will consider the single ratio B{fi —>■ 67) together with the leptogenesis 
constraint in Section [5l 
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Figure 2: Impact of iterative vs simplified procedure. Left plot: Simplified result for the 
ratio of branching ratios B{fi — > e'-f)/B{T — > fi'y), normalized to the one obtained with 
the iterative procedure. Right plot: Similarly for A12. 



It is also interesting to compare our elaborate iterative procedure of matching high- and 
low-energy parameters to a simpler procedure where we simply impose the weak-scale 
PMNS and neutrino mass parameters at the scale Aqut (Fig. [2] (left), corresponding to 
the MSSM with tan/3 = 2). It turns out that both procedures agree well for small scales 
M^. (This agreement is slightly worse for tan/3 = 10.) For large values > 10^^ GeV, 
deviations up to a few orders of magnitude can occur for some choices of parameters. 
It appears that this is usually due accidentally small branching ratios in one of the 
approaches. This is supported by the right plot in the Fig. [2], which shows a good 
agreement for the more fundamental flavour-violating quantity A12 up to the (expected) 
different overall normalization. 



4.4 CP asymmetries 

We are also in a position to illustrate and check numerically our qualitative discussion in 
Section [3] of the CP asymmetries relevant for leptogenesis. A thorough investigation of 
the baryon asymmetry follows in the next section. Fig. [3] shows the sum of the three CP 
asymmetries | Yli^il defined below ( l50l) . for the generic three- flavour case (left plot) and 
the CP asymmetry ei for the effective two- flavour case where only yi2 7^ (right plot). 
One can see clearly that in the latter case the dependence on is weak and slightly 
reciprocal. In fact this dependence is approximately proportional to In^ Aqut/^^ (black 
solid line) in agreement with expectations. The generic case is shown in the left plot for 
the SM (black solid) as well as the MSSM for tan (3 = 2 (red dot-dashed) and tan /3 = 10 
(blue dotted), with the remaining parameters given in the Figure caption. In contrast 
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Figure 3: Left plot: dependence of | the generic (3-flavour) case. Right 

plot: effective 2-flavour case. Normal hierarchy, 77jijgiit°st _ q q2 eV ; yi2 = 0.8, yi^ = 0.2, 
1/23 = 0.6 (3-flavour case), yi2 = 1 and 1/13 = 1/23 = (effective 2-flavour case). The PMNS 
phases have been taken to he 6 = a = (3 = vr/lO. Right plot: Effective two-flavour case; 
only ei is shown, on a linear scale. 



to the two-flavour case, there is strong dependence on M^, for > 10^^ GeV, when the 
contribution due to Y^, alone starts to dominate the RGEs (!25l).(!26ll. The precise form of 
the Ml, dependence is quite sensitive to the "angles" yij, but the roughly linear growth of 
I ej| in the regime of large Mj, appears to be general. However, the figure also clearly 
shows a strong dependence on the MSSM parameter tan/3 particularly for small M^. 
Indeed already for relatively small tan (3 = 10 the CP asymmetries can be more than an 
order of magnitude larger than in the SM. Moreover, in the case of the MSSM we observe 
a sign change at some scale ~ lO^^ q^v, which can be traced to the different relative 
signs between the terms on the right-hand sides of fl25|) and fl26l) . This example clearly 
demonstrates a rather dramatic dependence on details of the model. Finally, as in the 
case of the double ratios above, we investigated the impact of the iterative procedure 
compared to the simplified approach and found it to be generically small. Hence we feel 
justified to use the simplified procedure in Section [5] in order to save computer time. 



5 Leptogenesis in the extended MLFV Framework 
5.1 Preliminaries 

One of the most plausible mechanisms for creating the observed matter-antimatter asym- 
metry in the universe is leptogenesis, where a CP asymmetry generated through the 
out-of-equilibrium L- violating decays of the heavy Majorana neutrinos leads to a lepton 
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asymmetry which is subsequently transformed into a baryon asymmetry by {B + L)- 
violating sphaleron processes [H [9l [35] . 

Unfortunately, even in its simplest realization through the well-known seesaw mecha- 
nism [11], the theory has too many parameters. Indeed, as recalled in Section 2.4 in the 
framework of the standard model (SM) extended with three heavy Majorana neutrinos 
Ni{i = 1,2,3), the high-energy neutrino sector, characterized by the Dirac neutrino 
{m£)) and the heavy Majorana neutrino (M^j) mass matrices, has eighteen parameters. 
Of these, only nine combinations enter into the seesaw effective neutrino mass matrix 
M^^ m£) , thus making difficult to establish a direct link between leptogenesis and 
low-energy phenomenology [36] • Furthermore, there are six CP- violating phases which 
are physically relevant at high energies, while only three combinations of them are po- 
tentially observable at low energies. Therefore, no direct link between the sign of the 
baryon asymmetry and low-energy leptonic CP violation can be established, unless extra 
assumptions are introduced. 

Furthermore, additional assumptions are usually required to completely determine the 
high-energy neutrino sector from low-energy observables. Typical examples are the in- 
troduction of texture zeros in the Yukawa matrices or the imposition of symmetries to 
constrain their structure [37]. On the other hand, the heavy Majorana neutrino masses 
can range from the TeV region to the GUT scale, and the spectrum can be hierarchical, 
quasi-degenerate or even exactly degenerate [38]. Despite this arbitrariness, the heavy 
Majorana neutrino mass scale turns out to be crucial for a successful implementation 
of the leptogenesis mechanism. In particular, the standard thermal leptogenesis sce- 
nario with hierarchical heavy Majorana neutrino masses (Mi -C M2 < M3) requires 
Ml > 4 X 10^ GeV [39], if A'^i is in thermal equilibrium before it decays, or the more 
restrictive lower bound Mi > 2 x 10^ GeV [IQ] for a zero initial A'^i abundance. Since this 
bound also determines the lowest reheating temperature allowed after inflation, it could 
be problematic in supersymmetric theories due to the overproduction of light particles 
like the gravitino [4T] . 

It should be emphasized, that the above bounds are model dependent in the sense that 
they can be avoided, if the heavy Majorana neutrino spectrum is no longer hierarchical. 
For example, if at least two of the Ni are quasi-degenerate in mass, i.e. Mi ~ M2 , 
then the leptonic CP asymmetry relevant for leptogenesis exhibits the resonant behavior 
El ~ Mi/(M2 — Ml) [20], El]. In this case, it is possible to show that the upper bound on 
the CP asymmetry is independent of the light neutrino masses and successful leptogenesis 
simply requires Mi 2 to be above the electroweak scale for the sphaleron interactions to 
be effective. The quasi-degeneracy may also be achieved in soft leptogenesis where a 



26 



small splitting is induced by the soft supersymmetry breaking terms |42j . 

Another possibility which has been recently explored [T6l [T7] relies on the fact that ra- 
diative effects, induced by the renormalization group (RG) running from high to low 
energies, can naturally lead to a sufficiently small neutrino mass splitting at the lepto- 
genesis scale. In the latter case, sufficiently large CP asymmetries are generated. 

In the minimal seesaw scenario with only two heavy neutrinos the resulting baryon 
asymmetry in the SM turns out to be below the observed value [16]. On the other 
hand, this mechanism can be successfully implemented in its minimal supersymmetric 
extension (MSSM) 

It has been shown [18] that the above problems in the SM can be overcome in a more 
realistic scenario where the effects of a third heavy neutrino are also taken into account. 
In [18], leptogenesis was studied in the framework of a model where there are three right- 
handed neutrinos, with masses Mi ^ M2 -C M3. We will discuss this scenario below as 
a special limit of the MLFV framework. 

In view of the above, it is important to analyze leptogenesis in the extended MLFV 
framework, where CP violation is allowed both at high and low energies. In the MLFV 
scenario, right-handed neutrinos are assumed to be exactly degenerate at a high energy 
scale. In the limit of exact degeneracy, no lepton-asymmetries can be generated. How- 
ever, as previously emphasized, even if exact degeneracy is assumed at a high energy 
scale, renormalization group effects lead to a splitting of right-handed neutrino masses 
at the scale of leptogenesis, thus offering the possibility of viable leptogenesis in the 
extended MLFV framework. 



5.2 BAU in the RRL and Flavour Effects 

In leptogenesis scenarios the baryon asymmetry of the universe tjb arises due to non- 
perturbative sphaleron interactions that turn a lepton asymmetry into a baryon asym- 
metry. The predicted value of tjb has to match the results of WMAP and the BBN 
analysis for the primordial deuterium abundance [13] 

^^ = !^ = (6.3 ± 0.3) X 10-^°. (49) 
The lepton asymmetry is generated by out-of-equilibrium decays of heavy right-handed 
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Majorana neutrinos Ni and is proportional to the CP asymmetry e\ with 

j_ r{N,^Li<p)-r{Ni^Li^) 
El [r(iv, 0) + r(iv, ^ U 0)] ' 



(50) 



and / denoting the lepton flavour, that arises at one-loop order due to the interference 
of the tree level amplitude with vertex and self-energy corrections. 
A characteristic of the MLFV framework is that only admissible BAU with the help of 
leptogenesis is radiative and thereby resonant leptogenesis. The mass splittings of the 
right-handed neutrinos induced by the RGE are of similar size AM ~ 0{Mi,Yi^Yj) as 
the decay widths F ~ 0{Mi,Y,yYj). This is the condition of resonant leptogenesis. The 
CP asymmetry is for the lepton flavour / given by 

e\ = J2 '^((Y'^YMYMYJh) 9{Ml M], r|) (51) 

where g{Mf , jT"^) is an abbreviation for the full result given in [2T]. The total CP 
asymmetries Ei are obtained summing over the lepton flavours /. 

The baryon to photon number ratio rjB can then be calculated solving the Boltzmann 
equations for the lepton asymmetry and converting it into r]B using suitable dilution 
and sphaleron conversion factors. Which Boltzmann equation to use depends on the 
temperature scale at which leptogenesis takes place. We will follow a simplistic approach 
ignoring all subtleties generically coming into play in the intermediate regime between 
different mechanisms at work. Our main conclusions, however, will not be affected by 
this omission. We will simply divide the temperature scale into a region up to which all 
three lepton flavours have to be taken into account and a region above which the single 
flavour approximation works. 

Below some temperatur^ T^^^ ~ 10^ GeV [HI [45], |46], muon and tau charged lepton 
Yukawa interactions are much faster than the expansion H rendering the /i and r Yukawa 
couplings in equilibrium. The correct treatment in this regime requires the solution 
of lepton flavour specific Boltzmann equations. In the strong washout regime ?7_b is 
independent of the initial abundances and an estimate including flavour effects is given 

by m 

VB - -10-^ E e-^^^-^^^y^^^ e[ (52) 

j=l l=e,fj,,T 



^ We will chose T^^ ~ 10^" GcV in our analysis as an effective boundary between the unflavoured 
and 'fully flavoured' regimes, where we (respectively) neglect flavour and distinguish all three flavours. 
The main conclusions are, however, not affected by the precise choice. 
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with 

~ H{T = K'U) ^^^^ 

K,= Kl K^ = Y1 H{T = M,) ^ 17^ (54) 

l=e,fj,,T i=l 

where Mpi = 1.22 x 10^^ GeV and is the washout factor due to the inverse decay of 
the Majorana neutrino iVj into the lepton flavour /. The impact of lepton flavour effects 
on ?7b is discussed in [HI dSl SSI SSI SZ]- As we shall also see below, the inclusion of 
flavour effects generally leads to an enhancement of the resulting rjB- This is due to two 
effects: (1) the washout gets reduced since the interaction with the Higgs is with the 
flavour eigenstates only and (2) an additional source of CP violation arises due to lepton 
flavour speciflc CP asymmetries. 

For higher values of T > 10^^ GeV the charged lepton Yukawa couplings do not break the 
coherent evolution of the lepton doublets produced in heavy neutrino decays anymore. 
In this regime flavour effects can be ignored and an order of magnitude estimate is given 
by 

i=l l=e,fj,,T 

with K = ^^Ki. This agrees with a recent analytical estimate by [IE] up to factors of 
0{1) for the region of interest in parameter space, where the estimate of [IE] generally 
leads to a smaller efficiency and smaller rjB- We have also compared the analytical 



estimate of [H] and (I55 
the LeptoGen code 



with the numerical solution of the Boltzmann equations using 
For the relevant ranges of the input parameters the analytical 
estimate of (135]) and the full numerical solution agree quite well, whereas the estimate 
of [H] leads to an efficiency and rjB generally smaller by a factor of 5 to 10. This is 
shown in Fig. [HI These estimates, however, do not take into account the potentially 
large lepton flavour effects included in ( l52l) . 

Let us remark in passing that in the flavour independent region we are always in the 
strong washout regime, since 

K = K^ + K2 + K^ = —ti (Rd,R^) > ^ ^^-"^ ~ 50, (56) 

m* m* 

where m* = O{10~^). This inequality holds since the trace is linear function of the 
neutrino masses with positive coefficients, which reaches its minimum for yij = 0. We 
also made sure that the estimate (l52l) including flavour effects is applicable [19] and 



' http: / /www. ippp.dur.ac.uk/ ~teju/leptogen/ 
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checked that the inequahty 

Kl > 1 (57) 



is always satisfied for the points considered in the plots. Since both (l56l) and ( l57|) are 
satisfied, a simple decay-plus-inverse decay picture is a good description and the esti- 
mates (152|) and (!55l) independent of the initial abundances give a good approximation of 
the numerical solution of the full Boltzmann equations. 

We have performed the leptogenesis analysis specifically for the SM. We do not expect 
large deviations in the MSSM from the SM if the same Y^{M,y) and Ml{My) are given. 
The main differences come (1) from the CP-asymmetries, which now include contribu- 
tions from the supersymmetric particles, (2) from the washout, and (3) from conversion 
and dilution factors. The supersymmetric CP asymmetries have the same flavour struc- 
ture as in the SM and using [30] one can show that e^^^^^ ~ 2 e'^*^ for quasi-degenerate 
heavy neutrinos. We also expect the correction by the decay widths to be similar in size. 
Next, the washout in the strong washout regime is about a factor of -\/2 larger [l9] in the 
MSSM, whereas the dilution and sphaleron conversion factors stay almost unchanged. 
Concluding, we flnd that in the scenario considered the predicted values roughly satisfy 
^MSSM ^ i.5,^SM ^^Y^e same set of input parameters Yi,{My) and Mi{M^). The RGE 
induced values of Yy^M^) and Ml{My), however, are model dependent and lead to in 
general different Yy{My) and Ml{My) for the same boundary conditions at the GUT and 
low-energy scale, as discussed in Section [3^ Especially sensitive is the region My < 10^^ 
GeV where the CP asymmetries are dominantly generated by the tau Yukawa coupling, 
which is enhanced by a factor of tan (3 in the MSSM. Note also that in the MSSM, Tg'^ 
and TJq should be rescaled by a factor (1 + tan^ [3) to take account of the larger Yukawa 
couplings [59], which should make flavour effects even more prominent, § 



5.3 Two flavour limit 



As a flrst step we discuss the special case of yi2 being non-vanishing at the GUT scale 
and all other Uij = 0. This corresponds approximately to one of the scenarios considered 
in a recent study of radiative leptogenesis [18] with two right-handed neutrinos quasi- 
degenerate and a third right-handed neutrino decoupled Mi ~ M2 <^ M3. If only yi2 7^ 
the calculation of tjb proceeds in the same way since to a good approximation only z/^ and 
u]^ contribute to the CP asymmetry. The only difference comes from the enhanced wash- 
out. Since the third heavy neutrino is now also contributing, the lower bound on the 
washout K in (!56l) is in our case relatively enhanced by a factor (Am^^j^)^/^/ (Am^^^Y^'^ ~ 

^We thank S. Antusch for drawing our attention to this point. 
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Figure 4: Resulting rjs for the case in which only yi2 7^ (effective two flavour case) 
as a function of M^, for the normal hierarchy of light neutrinos: the orange crosses and 
red triangles show the unphysical limit setting the charged lepton Yukawas yg = in 
the renormalization group evolution with and without including lepton flavour effects in 
the calculation of 775, respectively. Setting the charged lepton Yukawas to their physical 
values, the blue circles and the green squares correspond to including and ignoring lepton 
flavour effects in the calculation, respectively. 



4 — 5. We have checked this correspondence for rjB also numerically. Ignoring flavour 
subtleties in leptogenesis for a moment, the CP violating effects due to renormalization 



group effects are induced only by the charged lepton yukawa couplings, see Section 13. 4[ 
and the total CP asymmetries for each heavy Majorana neutrino take the form [18] 



^1,2 



^1,2 



1 + D2,l 



(58) 



and 



tj — 



3|/2 lm{H2i) Re (F, 



u)13 



32 TT Hjj{H22-Hu] 



?>yl m.j{m.i + 777,2) ^/rr^lrl^ sinh(2?/i2) Re {U*2 Uri) 
64 TT {nil — ra2){m?, cosh^ yi2 + mim2 sinh^ 1/12) 



4 {H22 - Hiif In^ (M./Mout) 



TT 7T7^ COsh^ yi2 + 7772777,1 sluh^ yi2 



2 mj{m2 - nil) In {M^/Mi 



GUT J 



(59) 
(60) 



where Dj are regularization factors coming from the heavy Majorana decay widths. We 
immediately see that the total CP asymmetries only bare a very mild dependence on 
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Figure 5: rjB for the case in which only yi2 7^ (effective two ffavour function 
of yi2 (left) and mj,i (right) for the normal hierarchy. The black circles are obtained 
including lepton flavour effects and the red crosses are calculated ignoring them. 

the heavy Majorana scale. The almost negligible dependence on M^, has to be compared 
with the power-suppression in M^, in the hierarchical case (Mi -C M2 < M3). We find 
this expectation confirmed in Fig. HI where the resulting rjB is shown as a function of 
M,. 

Fig. m also nicely illustrates the relative importance of flavour effects in leptogenesis. If 
no cancellations occur, we find, that flavour effects generate an 77^ which is of the same 
order of magnitude (blue circles) , however almost always larger than the one calculated 
ignoring flavour effects (green squares). 

If we now consider the unphysical limit of setting = in the renormalization group 
running only, we find that the total CP asymmetries and tjb should vanish since no CP 
violation effects are induced by the RGE, see Section 13.41 We confirm this behavior in 
Fig. m (red triangles). A very different picture emerges once we include flavour effects. 
The relevant quantity for leptogenesis is then '^{{Y^Yj)ij(Y^)ii{Yj)ij) with no summation 
over the charged lepton index /. Although no total CP asymmetries are generated via 
the RG evolution in the limit Yg = 0, the CP asymmetries for a specific lepton flavour 
are no n- vanishing. Additionally, the resulting rjs now shows a M^, dependence which 
stems from the RGE contributions due to Y,y only, which are absent in the total CP 
asymmetries in the two flavour limit (orange crosses). 

All plots have been generated assuming a normal hierarchy of the light neutrino masses. 
We have checked that the results for the inverted hierarchy are similar, although rjs turns 
out to be generally smaller and below the observed value, in accordance with the findings 
of [IE]. Including flavour effects it is however still possible to generate a t]b of the correct 
order of magnitude. In Fig. [5] we additionally show the dependence of rjB on yi2 and 
m^i. We find that flavour effects enlarge the yi2 range where successful baryogenesis is 
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possible and slightly soften the upper bound on the light neutrino mass scale. The left 
panel even demonstrates that leptogenesis in the MLFV scenario is possible for a real R 
matrix. Then lepton flavour effects are essential for a successful leptogenesis |l5lll6]. 

5.4 General case 

Now we consider the general case with all three phases Uij non-vanishing. We have var- 
ied the parameters as described in Section HI The regularization of the resonant CP 
asymmetry by the Di turns out to be important for values of > 10~®, see Fig. |6l As 
seen there, in the regime where flavour effects are important we find an upper bound on 
the light neutrino mass of m^i < 0.2 eV in order to generate the right amount of 77^. 
Beyond the temperature scale where flavour effects play a role, no relevant bound can 
be found. This is due to the enhancement of the CP asymmetry which approximately 
grow linearly for values of > 10^^ GeV, see the discussion in Section 14.41 




Figure 6: (left) rjB for the general case with 0.01 < \yij\ < 1 as a function of rriyi 
(right). The hlack circles are obtained including lepton flavour effects whereas the red 
crosses are calculated ignoring them. The flavour blind results {red crosses) reach higher 
values due to the CP asymmetries growing as My gets bigger in this regime, (right) The 
total CP asymmetry |ei| for the general case with 0.01 < \yij\ < 0.8 as a function of 
Utot = (^12 + 2/13 + yis)^ input values that result in the right oder of magnitude of 
rjB- The red circles are obtained using the uncorrected CP asymmetries and the hlack 
squares include the corrections by the decay widths 

In Fig. [8], we compare different calculations of tjb'- 

• the flavour independent estimate of [H] used in Cirigliano et al. [13] (red boxes). 
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Figure 7: tjb for the general case with 0.01 < \yij\ < 1 as a function of M^,. The circles 
are obtained including lepton flavour effects and the red crosses are calculated ignoring 
them. 



• the numerical solution of the flavour independent Boltzmann equations using the 
LeptoGen package (black circles), 

• the recent estimate by Blanchet and Di Bari [4^ that includes flavour effects (green 
triangles), 

• the approximate expression of [19] given in (152!) that also includes flavour effects 
(brown crosses). 

We find that 

• the flavour blind estimate of [IB] used in Cirigliano et. al. [13] lies consistently 
below the numerical solution of the flavour independent Boltzmann equations. For 
Ml, > 10^^ GeV this turns out to be unimportant as flavour effects in this region 
are small and we confirm the increase oi rjB with M^, in this region found by these 
authors. 

• Potentially large flavour effects that have been left out in [13] generally enhance 
the predicted tjb, in particular for M^, < 10^^ GeV, in accordance with the existing 
literature. 

• Both flavour estimates and the numerical solution of the flavour independent Boltz- 
mann equations show solutions with ?7b of the measured order of magnitude without 
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imposing a stringent lower bound on the value of Mi 



The last finding is in contrast to the analysis of Cirigliano et. al. [13] which using the 
fiavour independent estimate of [H] finds a lower bound on My of (9(10^^GeV) as clearly 
represented by the red boxes in Fig. [HI The same qualitative conclusion holds for rjE 
using the RGE induced CP asymmetries in the MSSM. The tan (3 enhancement of the 
CP asymmetries as discussed in Section even facilitates the generation of an rjB of 
the right size. 

Our analysis that includes fiavour effects demonstrates that baryogenesis through lep- 
togenesis in the framework of MLFV is a stable mechanism and allows a successful 
generation of 77^ over a wide range of parameters. The absence of a lower bound on M^ 
found here has of course an impact on the LFV processes, which we will discuss next. 
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Figure 8: Different determinations of rjB for the general case with 0.01 < \yij\ < 1 
as a function of My. The black circles are obtained numerically solving the fiavour 
independent Boltzmann equations using the LeptoGen package, the green triangles and 
the brown crosses show estimates including fiavour effects of [17] and flS2]) . respectively. 
Finally the red boxes show the estimate of [IH] used in Cirigliano et. al [I^ which ignores 
flavour effects. 



In Fig. [9] we show B{fi ^ cy) vs. My for the parameter ranges described above and 
a lepton flavor violation scale of 1 TeV. We highlighted the points where successful 
baryogenesis is possible (black squares). We flnd that B{fi cy) can be made small 
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Figure 9: B{fi ^ 67) as a function of M^, for Alfv = 1 TeV. The black squares show 
points where a baryon asymmetry in the range 2 ■ 10^^° < 77^ < 10 ■ 10^^*^ is possible. 

enough to evade bounds from current and future experiments and one can have successful 
baryogenesis through leptogenesis at the same time. This is another finding of our paper 
which is in contrast to a recent analysis [13]. We will summarize the differences to [13] 
in the next paragraph. 



5.5 Comparison with [ jl3l 



Recently in an independent analysis Cirigliano, Isidori and Porretti [13] generalized 
MLFV formulation in [10] to include CP violation at low and high energy. Similarly 
to us they found it convenient to use for Yy the parametrization of Casas and Ibarra. 
They have also pointed out that in the MFLV framework the most natural is the resonant 
leptogenesis. 

On the other hand, these authors neglected flavour dependent effects in the evaluation of 
rjB, that we find in agreement with other authors to be important [HI [IHl US, HE], HTj . This 
has important consequences already at the qualitative level. Their qualitative discussion 
of the splittings of the at the see-saw scale is similar to ours and we agree with the 
main physical points made by these authors in this context. On the other hand, while 
we have demonstrated explicitely by means of a renormalization group analysis that a 
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successful RRL can be achieved, Cirigliano et al confined their analysis to parametrizing 
possible radiative effects in terms of a few parameters. In this context a new point 
made by us (see discussion Section 13.31) is that the coefficients q in (l20l) are in fact not 
independent of each other. Indeed the leading logarithmic contribution to q are related 
by the renormalization group. This can in principle increase the predictivity of MLFV. 
The three most interesting messages of [13] are 

• A successful resonant leptogenesis within the MLFV framework implies a lower 
bound > 10^2 GeV, 

• With Alfv = 0{1 TeV), this lower bound implies the rate for fi e'-f close to the 
present exclusion limit, 

• MLFV implies a specific pattern of charged LFV rates: -B(/i — » 67) < B{t ^ fi'j). 

For Ml, > 10^^ GeV, in spite of some differences in the numerics as discussed above, 
we basically confirm these findings. Unfortunately, for lower values of Mj, our results 
differ from theirs. In particular, as we have demonstrated in Fig. [8|, the observed value 
of T]B can be obtained for by several orders of magnitude below the bound in |13j . 
in accordance with other analyses of leptogenesis. Once is allowed to be far below 
10^^ GeV, Alfv = TeV) does not imply necessarily B{fi 67) close to the inclusion 
limit as clearly seen in Fig. [91 

One of the reasons for the discrepancy between our result with regard to Mj, and the 
one of [13] is the neglect of flavour effects in leptogenesis in the latter paper. Fig. M 
illustrates that flavour effects in leptogenesis matter. 

Concerning B{fi ^ 67) < B[t ^ H'j), we conflrm the result of [l^ in the limit of very 
small yi2, but as shown in Fig. [H this is not true in general, as also found in [26], [29] . 
Consequently, this hierarchy of charged LFV rates cannot be used as model independent 
signature of MLFV. 

6 Summary and Conclusions 

In this paper we have generalized the proposal of minimal flavour violation in the lepton 
sector of [10] to include CP Violation at low and high energy. While the deflnition 
proposed in [10] could be considered to be truly minimal, it appears to us too restrictive 
and not as general as the one in the quark sector (MFV) in which CP violation at low 
energy is automatically included [1] and in fact all flavour violating effects proceeding 
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through SM Yukawa couphngs are taken into account [2|. The new aspect of MLFV in 
the presence of right-handed neutrinos, when compared with MFV, is that the driving 
source of flavour violation, the neutrino Yukawa matrix Y^, depends generally also on 
physics at very high scales. This means also on CP violating sources relevant for the 
generation of baryon-antibaryon asymmetry with the help of leptogenesis. The first 
discussion of CP violation at low and high energy has been presented in [13]. Our 
conclusions for M^, > 10^^ agree basically with these authors. However, they differ in an 
essential manner for lower values of M,^. 

The main points of our paper have been already summarized in the introduction. There- 
fore it suffices to conclude our paper with the following messages: 

• A new aspect of our paper is the realization that in the context of MLFV the only 
admissible BAU with the help of leptogenesis is the one through radiative resonant 
leptogenesis (RRL). Similar observations have been made in [13]. In this context 
our analysis benefited from the ones in [16l [TTJ [20], El]. The numerous analyses 
of leptogenesis with hierarchical right-handed neutrinos present in the literature 
are therefore outside the MLFV framework and the differences between the results 
presented here and the ones found in the literature for Mi <^ M2 <^ M3 can 
be used to distinguish MLFV from these analyses that could be affected by new 
flavour violating interactions responsible for hierarchical right-handed neutrinos. 

• We have demonstrated explicitely within the SM and the MSSM at low tan /3 that 
within a general MLFV scenario the right size of 77^ can indeed be obtained by 
means of RRL. The important property of this type of leptogenesis is the very 
weak sensitivity of 77^ to the see-saw scale so that for scales as low as 10^ GeV 
but also as high as 10^'^ GeV, the observed tjb can be found. 

• Flavour effects, as addressed by several authors recently in the literature [HI [191 
H5l [IHl [17], play an important role for < 10^° GeV as they generally enhance 
rjE- Moreover, they allow for a successful leptogenesis within MLFV even when 
the i?-matrix is real (left panel of Fig. [5]). 

• As charged LFV processes, like /i — > 67 are sensitive functions of M^, while rjB 
is not in the RRL scenario considered here, strong correlations between the rates 
for these processes and 77^, found in new physics scenarios with other types of 
leptogenesis can be avoided. 

• Except for this important message, several of the observations made by us with 
regard to the dependence of charged LFV processes on the complex phases in the 
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matrix R and the Majorana phases have been aheady made by other authors in 
the rich hterature on LFV and leptogenesis. But most of these analyses were done 
in the context of supersymmetry. Here we would like to emphasize that various 
effects and several patterns identified there are valid also beyond low energy super- 
symmetry, even if supersymmetry allows a definite realization of MLFV provided 
right-handed neutrinos are degenerate in mass at the GUT scale. 

• One of the important consequences of the messages above is the realization that 
the relations between the flavour violating processes in the charged lepton sector, 
the low energy parameters in the neutrino sector, the LHC physics and the size 
of ?7b are much richer in a general MLFV framework than suggested by [TOl [13]. 
Without a specific MLFV model no general clear cut conclusions about the scale 
Alfv on the basis of a future observation or non-observation of /i — > 67 with the 
rate O{10~^^) can be made in this framework. 

• On the other hand we fully agree with the point made in that the observation 
of /i — s> 67 with the rate at the level of 10"^'^, is much easier to obtain within 
the MLFV scenario if the scales Alfv and M^, are sufficiently separated from each 
other. We want only to add that the necessary size of this separation is sensitive 
to the physics between Mz and Agut, Majorana phases and CP violation at high 
energy. In this manner the lepton flavour violating processes, even in the MLFV 
framework, probe scales well above the scales attainable at LHC, which is not 
necessarily the case within MFV in the quark sector. 

• Finally, but very importantly, MLFV being very sensitive to new physics at high 
energy scales, does not generally solve possible CP and flavour problems. This 
should be contrasted with the MFV in the quark sector, where the sensitivity to 
new physics at scales larger than 1 TeV is suppressed by the GIM mechanism. 



Note 

During the preparation of this revised version, one of us (S.U.) has investigated para- 
metric dependences in the present scenario for the case of a real R in more detail [60] . 
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A Iterative solution of the renormalization group 
equations 

The goal of our numerical analysis of Section 4 is to determine the neutrino Yukawa 
matrix Y^, and the masses of the right-handed neutrinos at the scale Mj, taking into 
account the constraints on the masses and mixings of light neutrinos measured at low 
energies and imposing the GUT condition characteristic for the MLFV 

Mi(Agut) = M2(Agut) = M3(Agut). (61) 
As discussed in Section 2 the latter condition implies 

Re(i?(AGUT)) = 0, (62) 
but Im(i?(AGUT)) niust be kept non-zero in order to have CP- violation at high energy. 

The RG evolution from Agut down to generates small sphttings between Mi{M^) and 
a non- vanishing Re(i?(M^)), both required for the leptogenesis. As the splittings between 
Mi{My) turn out to be small, we integrate the right- handed neutrinos simultaneously 
at /i = My imposing, up to their splittings, 

Mi(M,) ^ M2(M,) ^ M^{M,) ^ M,. (63) 

In view of various correlations and mixing under RG between different variables we reach 
the goal outlined above by means of the following recursive procedure: 

Step 1 

We associate the values for the solar and atmospheric neutrino oscillation parameters 
given in ( H0l) - (HT1) with the scale /i = Mz and set 6*13 and the smallest neutrino mass 
^lightest particular values corresponding to ^ = Mz- 
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step 2 

For a chosen value of M^, the RG equations, specific to a given MLFV model, are used 
to find the values of the parameters of Step 1 at /i = M^. For instance we find m^(Mjy) 
and similarly for other parameters. 

Step 3 

We choose a value for Agut and set first m^(AGUT) = m^{Mi,) and similarly for other 
parameters evaluated in Step 2. Setting next 

Ml (Agut) = M2(Agut) = M3(Agut) = (64) 

and choosing the matrix R, that satisfies (j62l) . allows also to construct Yi,{Ag\jt) by 
means of the parametrization in (ITTll . 

Step 4 

Having determined the initial conditions for all the parameters at /i = Agut we use 
the full set of the RG equations [53] to evaluate these parameters at M^. In the range 
Ml, < fi < Agut we use 

niM = -v'Yj{fi)M-\fi)YM. (65) 

The RG effects between Agut and M^, will generally shift m^{M^) to new values 

<(M,) = m';{M,) + Am^ (66) 

with similar shifts in other low energy parameters. If these shifts are very small our goal 
is achieved and the resulting Y^{Mi,) and Mi{My) can be used for lepton flavour violating 
processes and leptogenesis. If the shifts in question are significant we go to Step 5. 

Step 5 

The initial conditions at /i = Agut are adjusted in order to obtain the correct values for 
low energy parameters at = My as obtained in Step 2. In particular we set 

mr(AGUT) = mr(M,) - Am^ (67) 

with AmJ' defined in ( l66l) . Similar shifts are made for other parameters. If the condition 
fl55]) is not satisfied in Step 4, the corresponding shift in (1641) should be made. Choosing 
R as in Step 3 allows to construct an improved yj,(AGUT)- Performing RG evolution with 
new input from Agut to M^ we find new values for the low energy parameters at My 
that should now be closer to the values found in Step 2 than it was the case in Step 4. 
If necessary, new iterations of this procedure can be performed until the values of Step 
2 are reached. The resulting Yy{My) and Mi[My) are the ones we were looking for. 
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B Basic Formulae for 



B.l Preliminaries 

In what follows we will present two generalizations of the formulae for Ajj in [10] in 
the approximation of degenerate right-handed neutrinos. We have checked that the 
splitting of by RGE has very small impact on these formulae. All the expressions 
for Aij are meant to be valid at M,^. Similar formulae have been given for instance in 
[231 [261 [271 [25l [3T] . but we think that the formulae given below are more transparent. 

In order to obtain transparent expressions for Aij it is useful to introduce the mass 
differences 

S21 = m^2 - m^i, S31 = m^s - m^i, (68) 
'^'21(^12) = ^21 cosh(2yi2), 531(^12) = + m^i{l - cosh(2?/i2)) (69) 
and collect the dependence on Majorana phases in the following two functions 

Fi(a, (3) = e-^("-^) + 2ic^ sin(a - /3), (70) 

^2(0, P, S) = sisc^ cos(a — P + S) + ics sin(a — P) + si^s^ cos(a — P — S). (71) 



B.2 R real and Ui^ complex 

In this case we find 

A^, = -^{s c + e-''sM, (72) 
M 

Are = ^(-^ C 521 + e-^'sM, (73) 

Ar, = ^{-c%, + S,,), (74) 
where we have neglected terms (^(sia), whenever it was justified. 

For the CP conserving cases 6 = 0, vr, these formulae reduce to the formulae (24) of [TO] 
which represent the case where both matrices, R and U^, are real. We note that in the 
presence of a real matrix R, the Aij do not depend on the Majorana phases and A,-^ 
does not depend on 6. 
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B.3 R and Ui^ complex 

Allowing for one additional phase yu in R, we find the generalization of (172l) -(!7 ^ that 
includes CP violation both at low and high energies represented by 5 7^ 0, tt and X12, yi2 7^ 
0, respectively 



A^e = (^s c 521(^12) + e~*'^si353i(2/i2) + i^Jmyimy2 sinh(2?/i2)Fi (a, , (75) 

A^e = (y-s c 521(2/12) + e~*'^si353i(yi2) - iy/m^im^2 sinh(2?/i2)Fi(Q;, /?)) , (76) 

A^^ = ^ (^-c^ ^21(2/12) + ^31(2/12) + 2iy/m^im^2 sinh(2yi2)F2(a, /3, 5) j . (77) 

For yi2 = ([75])-([77D reduce to ([72])-([74D. We note that relative to ([72D-([74D there is 
an additional dependence on the difference of Majorana phases a — /5, collected in the 
functions Fi and F2 that disappears for yi2 = 0. This means that for yi2 very close to 
zero Majorana phases in /j — *■ Ij'j decays do not matter but can be important already 
for small yi2. 

Indeed, the Aj^'s are very sensitive to y^ and the values of Aij can be enhanced by 
several orders of magnitude [SSlEllEHlEHlETlEalEIllal] relative to the case of i? = 1, 
even for yi2 = Indeed as seen in (175l) - (!77jl . the Aij depend exponentially on the yi2 

and moreover for yi2 7^ 0, they do not only depend on the neutrino mass differences but 
also on ^Jmvlmy2 which can be much larger than Arriij. Thus including a non- vanishing 
phase in R can have in principle a very strong impact on the analysis of [10] as also 
discussed in 



The large enhancement of A^- in the case of a complex R is analogous to the large 
enhancement of B{Bd,s /^^A^~) for large tan/3. In the latter case the presence of new 
scalar operators lifts the helicity suppression of the branching ratios in question. In the 
case of Aij the appearance of a new mass dependence mirrij in addition to rrii — rrij has 
a similar effect provided y/ffhirnj mi — rrij . 
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